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Shell evolution is studied in the neutron-rich silicon isotopes 36,38,40gj ygi^g neutron single-particle 
strengths deduced from one-neutron knockout reactions. Configurations involving neutron excita¬ 
tions across the A = 20 and A = 28 shell gaps are quantified experimentally in these rare isotopes. 

Comparisons with shell model calculations show that the tensor force, understood to drive the col¬ 
lective behavior in ^^Si with A = 28, is already important in determining the structure of ^°Si 
with A = 26. New data relating to cross-shell excitations provide the first quantitative support for 
repulsive contributions to the cross-shell T = 1 interaction arising from three-nucleon forces. 


The atomic nucleus is a fermionic many-body quan¬ 
tum system composed of strongly-interacting protons 
and neutrons. Large stabilizing energy gaps, separating 
clusters of single-particle states, provide the cornerstone 
for the nuclear shell model, one of the most powerful 
tools available for describing the structure of atomic nu¬ 
clei. In the simplest version of the shell model, empirical 
shell gaps at the magic nucleon numbers 2, 8 , 20, 28, 50, 
82, and 126 are reproduced when assuming that the nu¬ 
cleons experience, predominantly, a mean-field potential 
with an attractive one-body spin-orbit term. 

In rare isotopes, with imbalanced proton and neutron 
numbers, significant modifications have been observed. 
Here, new shell gaps develop and the conventional gaps 
at the magic numbers can collapse. Understanding this 
observed evolution is key to a comprehensive description 
of atomic nuclei across the nuclear chart. Detailed stud¬ 
ies of the evolution of shell structure with proton number 
{Z) or neutron number (A), e.g. [T], probe the effects 
of particular components of the complex interactions be¬ 
tween nucleons: such as the spin-isospin [2] and tensor [3] 
two-body terms, and three-body force terms [SIS]. The 
need to include such terms in the nuclear interaction has 
been demonstrated by their robust effects that become 
amplified at large isospin nig and, without which, fea¬ 
tures such as driplines and shell structure may not be 
reproduced. Clearly, a full treatment of the nuclear force 
from its underlying QCD degrees of freedom is very chal¬ 
lenging, and experimental data is essential in helping to 
identify the most important degrees of freedom responsi¬ 
ble for driving the evolution of nuclear properties. 

Here, we present data for the silicon {Z = 14) isotopic 
chain, a region of the nuclear chart were rapid shell evo¬ 
lution is at play. ^"*81 is known to exhibit closed-shell be¬ 


havior while ^^Si shows no indication of an A = 28 shell 
gap m- Hitherto, observations on the neutron-rich sil¬ 
icon isotopes, dominated by measurements of collective 
observables, have been reproduced by large-scale shell 
model calculations using phenomenological effective in¬ 
teractions 1 ^. To assess the theoretical description of 
the evolving shell structure it is also critical to investigate 
these nuclei using single-particle observables, such as the 
energies and single-particle (spectroscopic) strengths of 
states involving the active orbitals at shell gaps. 

This Rapid Communication reports a first experimen¬ 
tal investigation of observables that reflect single-neutron 
degrees of freedom in the 36,38,40gj isotopes. Extraction 
of the presented cross-sections from the data, collected 
in the measurement reported in Ref. m, required the 
development of novel analysis strategies. The results go 
beyond those of Ref. [10] , and are interpreted here within 
a common theoretical framework. Exclusive one-neutron 
knockout cross sections, measured using y-ray tagged 
neutron removal reactions from 36,38,40gj projectiles, are 
used to identify and quantify configurations that involve 
neutron excitations across the A = 20 and N — 28 shell 
gaps. Specifically, the partial cross sections to the lowest- 
lying 7/2“ and 3/2“ states, involving the diminishing 
A = 28 gap, and 1/2+ and 3/2+ states, involving the 
A = 20 shell gap, are measured and compared to cal¬ 
culations using shell-model spectroscopic strengths and 
eikonal reaction theory. The results (i) track the evolu¬ 
tion of the neutron fY /2 and P 3/2 orbitals at the A = 28 
shell gap, and (ii) quantify the little-explored neutron ex¬ 
citations, from the c? 3/2 and Si /2 sd-shell orbitals, across 
the A = 20 gap. 

The experiment was performed at the Coupled Cy¬ 
clotron Facility of the National Superconducting Cy- 
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FIG. 2. (color online) Momentum distribution of the ®®Si 
residues produced in the ground state and any isomers fol¬ 
lowing one neutron knockout from ^®Si. The curves show a 
fit using a linear combination of the calculated distributions 
for removal from the / 7/2 and d 3/2 orbits. 
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FIG. 1. Experimental level schemes for low-lying states in 
®®Si, ^^Si and ^®Si, compared with shell model calculations 
using the SDPF-MU effective interaction [lOj . 

clotron Laboratory at Michigan State University. Sec¬ 
ondary beams of 36.38,40gj^ produced by fast fragmenta¬ 
tion of a “^^Ca primary beam, impinged on a beryllium 
target with energies of 100, 95 and 85 MeV/u, respec¬ 
tively. The one-neutron knockout residues were detected 
and identified on an event-by-event basis. Prompt 7 - 
rays, emitted in-flight from de-excitation of the knock¬ 
out residues, were detected with the GRETINA array 
m surrounding the target position, and were Doppler- 
corrected event-by-event. 

The level schemes of the knockout residues were con¬ 
structed based on 77 coincidences, energy sums, and in¬ 
tensity balances. These are summarized in Fig. Spin- 
parity assignments were made with the aid of the paral¬ 
lel momentum distributions of the residues in compari¬ 
son with theoretical distributions calculated in an eikonal 
model according to the formalism of Ref. [T^] . Full de¬ 
tails of the experiment, data analysis and the spin-parity 
assignments can be found in Ref. uni- 

The knockout cross sections to negative parity states, 
i.e. removal from the neutron / 7/2 and P 3/2 orbitals in 
36,38,40gj^ map their spectroscopic strengths. The experi¬ 
mental and calculated cross sections are listed in Table HI 
Details of the shell model calculations used can be found 
in Ref. (TU]. Determining the partial cross sections is 
challenging in some cases. For example, population of 
the ^^Si(7/2)") ground state is hindered by the presence 
of a 3/2+ isomer, expected to be strongly populated but 
which cannot be tagged with prompt in-beam 7 spec¬ 
troscopy. We use instead the ^®Si residue momentum 
distribution to extract the population fraction. Fig. 


TABLE I. Experimental {cj^xp) and calculated {(Jth) one- 
neutron knockout cross sections to the lowest 7/2“ and 3/2- 
states in the mass Ares residues. The ath use the shell-model 
spectroscopic factors C^S and their center-of-mass correc¬ 
tion, {Aproj/Ares)^, and the calculated eikonal model single¬ 
particle cross sections asp. All cross sections are in millibarns. 


TTT 

^.f 

■^proj 

-^res 

sp 

SDPF-MU 
C^S ath 

SDPF-U 
C^S ath 

(^exp 

7/2- 

36 

35 

15.7 

1.71 

29.2 

1.73 

29.5 

23(6) 


38 

37 

15.2 

2.81 

46.3 

2.85 

46.6 

47(9) 


40 

39 

15.0 

3.19 

51.8 

3.33 

53.9 

49(7) 


42 

41 

15.0 

2.73 

46.6 

3.70 

64.1 


3/27 

36 

35 

17.8 

6.i3 

2.6 

0.09 

1.8 

8(3) 


38 

37 

20.1 

0.11 

2.4 

0.02 

0.4 

9(7) 


40 

39 

21.0 

0.90 

21.3 

0.51 

12.2 

29(20) 


42 

41 

26.0 

1.72 

49.4 

0.03 

1.0 


3/27 

38 

37 

18.8 

0.27 

5.5 

0.27 

5.6 

7(3) 


40 

39 

19.0 

0.08 

1.8 

0.17 

3.7 



42 

41 

22.0 

0.19 

4.6 

1.02 

30.5 



Experimental excitation energy not known, ccsp is calculated 
with the SDPF-MU shell model energy. 


shows the ^^Si parallel momentum distribution after the 
subtraction of all events that decay by prompt 7 emis¬ 
sion. Overlaid is a linear combination of the theoretical 
distributions for neutron removal from the / 7/2 and (^ 3/2 
orbits together with the resulting fit minimization, 
giving a / 7/2 fraction of 0.45(10). The ground state cross 
section is estimated from this fraction of the knockout 
reaction events with no prompt 7 decay. 

The ^''Si(3/2/) and ^®Si(7/2/) states are nanosecond 
isomers. As a result, their depopulating transitions have 
broad, asymmetric peak shapes due to the larger uncer¬ 
tainty in the position and velocity of the decaying frag¬ 
ment and corresponding degradation of the Doppler re¬ 
construction. This lifetime effect is incorporated into the 
GEANT4 [T3] simulation and a best-fit lifetime is ob- 
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tained with a maximum likelihood method. An example 
is shown in Fig. The proximity of these peaks to the 
7 -ray detection threshold results in a dependence of the 
extracted peak intensity on the assumed lifetime. This 
dependence is shown in the lower panel of Fig. Ufa) , and 
makes the major contribution to the uncertainty in these 
peak intensities. The effect of the lifetime on the peak 
shape depends on the polar angle of the emitted 7 ray. 
We can confirm that the simulation reproduces this de¬ 
pendence by dividing the array into three rings centered 
near 50, 65 and 90 degrees (labeled front, middle, and 
backward in Fig. |^, and comparing the fit in each ring. 
This comparison, in Fig. [^b), shows satisfactory agree¬ 
ment. 





Energy (keV) 


FIG. 3. (color online) (a) Maximum likelihood fit to the life¬ 
time of the state in ^®Si decaying by a 172 keV 7 ray. Only 
statistical uncertainties are shown. The lower panel shows the 
influence of the uncertainty in lifetime upon the uncertainty in 
the number of decays, where the error band shows the uncer¬ 
tainty in the ht. (b) The resulting spectra in the front, middle 
and backward rings of GRETINA using the best-ht lifetime. 
The shaded gray area indicates the continuous background 
from target breakup, the blue area indicates the Compton 
continuum of all higher-lying peaks, and the magenta curve 
indicates the total ht | 10 |. 

The large uncertainty for the ^®Si(3/2(") state in Ta¬ 
ble U is due to several observed transitions which were 
not placed in the level scheme, introducing ambiguity in 
the subtraction procedure described above. The quoted 
uncertainty includes the range of possible level schemes 
which are consistent with the data. Further, since the 
second 3/2“ state was not identified, the value shown 
provides only a lower limit on the bound P 3/2 strength. 
The stated ^^’^®Si(7/2“) cross sections assume that pop¬ 
ulation of the predicted 5/2/ states is small compared to 
other sources of uncertainty (the shell model strengths 
predict cross sections of order 1 mb). The measured 


FIG. 4. (color online) One-neutron knockout reaction partial 
cross sections to bound final states with J'" = 7/2“ and 3/2“, 
as a function of the neutron number of the silicon projectile. 
The arrow, for N = 26, indicates a lower limit. Theoreti¬ 
cal predictions using the SDPF-MU (solid lines) and SDPF- 
U (dashed lines) effective interactions are also shown. The 
dotted line (SDPF-MU-NT) results when the tensor compo¬ 
nent of the cross-shell sd-fp interaction of SDPF-MU is set to 
zero. All theoretical curves have been scaled by the observed 
knockout reaction systematics (see text). 


and theoretical cross sections (for the SDPF-MU [5] and 
SDPF-U [3 shell model effective interactions) are shown 
in Fig. 1^ In the region of ^^Si, the tensor component 
of the interaction has been proposed as an important 
driving force for shell evolution HIS], and so we investi¬ 
gate this with a third set of calculations—denoted SDPF- 
MU-NT—obtained by removing the tensor part of the 
cross-shell sd-fp interaction of SDPF-MU. All theoreti¬ 
cal cross sections are scaled by an empirical quenching 
factor R{AS) obtained from a fit to the knockout reac¬ 
tion systematics [iniiiii- 

The agreement between the measured 7/2“ state cross 
sections (shown in blue) and both the SDPF-MU and 
SDPF-U calculations is excellent. We see that the ef¬ 
fect of the tensor force, as discussed in [ 2 , becomes im¬ 
portant m already around '^°Si. In contrast, the 3/2 
state cross sections (shown in red) are markedly under¬ 
predicted. This finding is consistent with previous mea¬ 
surements using one-neutron knockout, from 30’32]y[g 
and ^^Mg |T7] . as well as a {t,p) transfer measurement 
populating states in ^^Mg [18]. In each of these cases, an 
excess of P 3/2 strength was seen, relative to shell model 
predictions, while the / 7/2 strength was generally consis¬ 
tent with the shell model. The fact that this discrepancy 
is observed for different reaction mechanisms, and only 
for a particular orbit, suggests that this is a structure 
effect and not related to any systematic defect of the re¬ 
action theory. As can be seen from the dotted line in 
Fig. a the tensor force does not appear to have much 
effect in ^®Si and ^®Si, and so the P 3/2 discrepancy likely 





















4 


TABLE II. Experimental {aexp) and calculated (ath) one- 
neutron knockout cross sections to the lowest 3/2+ and 1/2+ 
states in the mass Ares residues. The ath use the shell-model 
spectroscopic factors C^S and their center-of-mass correc¬ 
tion, (AprojlAres)^, and the calculated eikonal model single¬ 
particle cross sections asp. All cross sections are in millibarns. 






SDPF-MU 

SDPF-U 


JJ 

■^proj 

A 

■^res 

sp 

C'^S 

ath 

C^S 

ath 

^exp 

3/2+ 

36 

35 

14.2 

3.07 

46.0 

2.61 

39.2 

29(6) 


38 

37 

13.8 

2.79 

40.6 

2.19 

28.4 

19(2) 


40 

39 

13.3 

2.31 

32.4 

1.69 

19.9 

14(2) 


36 

35 

21.1 

0.96 

21.3 

1.00 

22.3 

13(1) 


38 

37 

20.7 

0.80 

17.5 

0.97 

19.5 

10 ( 1 ) 


40 

39 

22.7(3 

0.53 

12.6 

0.72 

14.1 



Experimental excitation energy not known, asp is calculated 
with the SDPF-MU shell model energy. 


has origins elsewhere. 

To clarify the N = 20 shell closure we also consider the 
removal of neutrons from the ^ 3/2 and si /2 sd-shell or¬ 
bitals, populating positive-parity final states. The cross 
sections for population of bound 3/2+ and 1/2+ states 
are listed in Table The large uncertainty for the 
^®Si(3/2+) state yield is due to the same isomer effect 
as was discussed for the 7/2]" state. The measured and 
theoretical (SDPF-MU and SDPF-U) cross sections are 
compared in Table |TT| Both model calculations, which in¬ 
clude Ip — Ih excitations from the sd-shell in the wave 
functions of the residual nuclei, over-predict the strength 
of transitions from these orbits. 

It is very likely that the theoretical over-prediction 
of sd-shell strength and the aforementioned under¬ 
prediction of /p-shell strength are related, reflecting 
unaccounted-for excitations across the = 20 shell 
gap in the ground states of the projectiles. Indeed, the 
present calculations for the projectile ground states were 
performed in a OHuj model space in which the neutron 
sd-shell orbits are fixed and fully occupied. So, it is ev¬ 
ident that the assumed occupation of sd-shell orbits is 
too high. In the Monte Carlo shell model calculations 
of Ref. |19] . that allowed an arbitrary number of neu¬ 
tron particle-hole excitations from the sd-shell into the 
lower fp shell, the results, in ^®’^®Si, were an average 
excess of approximately 0.3—0.4 neutrons compared to 
normal filling. This reduced sd strength (and additional 
fp strength) would bring the shell-model predictions into 
better agreement with the present data, with the excep¬ 
tion of the large 3/2+ strength of SDPF-MU. 

Finally, the newly-measured energies of the 3/2+ and 
1 / 2 + hole states provide guidance for shell-model ef¬ 
fective interactions that include excitations across the 
iV = 20 shell gap. Figure shows the experimental ener¬ 
gies of the 3/2+ (1/2+) states relative to the 7/2f states, 
indicative of the / 7/2 to ^ 3/2 (si/ 2 ) shell gap. For refer¬ 
ence, we also show the shell-model spectroscopic factors 


for populating these states by one-neutron removal. The 
experimental data indicate that both gaps shrink as neu¬ 
trons are added from TV = 19 to 25, while SDPF-MU 
predicts a flat trend and SDPF-U predicts an increase of 
these gaps. 

These qualitiatively different predictions can largely 
be attributed to a difference in the cross-shell neutron- 
neutron (T = 1 ) interaction. Figure shows selected 
monopole (i.e. angle-averaged) terms of the SDPF-U 
and SDPF-MU interactions. While both interactions 
have similar sd and fp monopoles and are successful in 
reproducing the spectroscopy of the region within the 
Ohuj model space, the more-attractive SDPF-U cross¬ 
shell monopoles over-bind the neutron sd orbits as neu¬ 
trons are added to the fp shell, leading to the observed 
trend. This discrepancy highlights a key difference be¬ 
tween the two interactions. In SDPF-U, due to insuf¬ 
ficient experimental data, the cross-shell part of the in¬ 
teraction was left as essentially the two-body G matrix. 
On the other hand, the cross-shell component of SDPF- 
MU was generated from the schematic potential Vmu [20] 
which allowed—by incorporating information from data 
closer to stability—the approximate inclusion of the re¬ 
pulsive contribution of three body forces to the effective 
T = I two body interaction. This same repulsive T = I 
effect has been shown to be robust consequence of the 
Fujita-Miyazawa process which is crucial in reproducing 
the oxygen dripline [5|. We note that a more recent ver¬ 
sion of SDPF-U |9], developed to allow neutron excita¬ 
tions across the = 20 gap, in fact produces significant 
improvement over the original SDPF-U energies |21j . 



FIG. 5. (color online) Energies of 3/2+ and 1/2+ intruder 
states in the silicon isotopes as a function of mass number. 
The points indicate experimental data, while the lines con¬ 
nect theoretical calculations. The vertical bars show the shell- 
model spectroscopic factors of the states. 

In conclusion, we have exploited one-neutron knockout 
reactions to probe the evolution of the / 7/2 and P 3/2 spec¬ 
troscopic strength in neutron-rich silicon isotopes. State- 
of-the-art shell-model interactions describe the trends of 
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FIG. 6. (color online) Comparison of selected T = 1 
monopoles terms of the SDPF-U and SDPF-MU effective in¬ 
teractions, evaluated at ^ = 42. 

the data but underestimate the role of the p ^/2 orbital. 
We confirm that the tensor force is necessary to describe 
the evolution of the f’j /2 strength, and show that it is 
already important aX N = 26. The observed excess of 
P 3/2 strength relative to shell-model predictions indicates 
that the N = 28 shell gap may be reduced even more 
than present calculations suggest. Neutron cross-shell 
excitations across the iV = 20 shell gap were identified 
and quantified for the first time from the observation 
of positive-parity final states. The shell-model interac¬ 
tions considered (SDPF-U and SDPF-MU) over-predict 
the measured ^ 3/2 and Si /2 neutron removal yields, point¬ 
ing to the deficiency of the applied model space trunca¬ 
tions. We have also identified the energies of neutron- 
hole states which depend strongly on previously uncon¬ 
strained neutron-neutron monopole interactions. A com¬ 
parison of shell-model predictions indicates the impor¬ 
tance of three-body forces in the evolution of structure 
in this region. 
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